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ABSTRACT 

The  relation  between  the  site  of  energy  deposition  and  the  site  of  its  biological  action 
is  an  important  question  in  radiobiology.  Even  at  77*K,  evidence  is  clear  that  these 
two  sites  must  be  separated  since  energy  deposition  is  random  hut  specific  products 
are  formed.  Several  processes  that  may  contribute  to  this  separation  are:  1)  hole 
migration  and  stabilization  through  deprotonation  to  give  neutral  oxidation  product 
radicals;  2)  electron  trapping  and  transfer  to  form  specific  radical  anions,  possibly 
followed  by  protonation  to  give  neutral  reduction  product  radicals;  and  3)  recombina¬ 
tion  of  spatially  separated  charges  or  radicals.  These  microscopic  processes  will  be 
reviewed  critically  in  an  analysis  using  electron  paramagnetic  resonance  spectroscopy 
(EPR)  evidence  for  and  against  long-range  transfer  of  energy  and/or  charge  in  frozen, 
hydrated  DNA. 

INTRODUCTION 

A  strongly  held  tenet  of  radiobiology  is  that  damage  of  DNA  mediates  many  of  the 
cellular  effects  of  ionizing  radiation.  Whether  die  damage  results  from  attack  on 
DNA  by  OH  (hydroxyl  free  radical)  or  other  radicals  formed  at  a  distance  from  a 
DNA  molecule  (indirect  effect)  or  from  energy  deposited  directly  in  a  DNA  molecule 
by  radiation  (direct  effect)  has  been  the  subject  of  debate  for  many  years  III. 
Growing  evidence  indicates  that  die  direct  effect,  defined  to  include  deposition  of 
energy  in  DNA  and  water  molecules  closely  associated  with  DNA  (die  hydratioo 
layer  or  layers)  plays  a  critically  important  role  in  cellular  responses  to  radiation 
/23/.  It  follows  that  investigations  of  die  localization  of  energy  and  charge  in 
irradiated  DNA  are  urgently  needed. 

Experimental  studies  have  shown  convincingly  that  the  radiation  chemistry  of  DNA 
is  critically  dependent  on  experimental  conditions.  The  type  of  salt  and  its  concentra¬ 
tion  141,  and  the  state  of  DNA  hydration  /5,6,7/  are  among  the  variables  that  effect 
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DNA  chemistry.  EPR  spectroscopy  has  provided  a  most  valuable  experimental  tool 
to  develop  insight  into  the  fundamental  processes  involved  in  DNA  damage  reactions. 
EPR  spectra  ranging  from  a  0.9mT  "singlet”  spectrum  to  a  14mT  octet  spectrum  were 
reported  for  irradiated  DNA  as  early  as  1972  /8A  Sorting  out  the  factors  contributing 
to  spectral  variation  has  taken  the  efforts  of  many  outstanding  investigators.  A  major 
contribution  in  elucidating  DNA  radiation  chemistry  at  low  temperature  was  the 
preparation  by  Rupprecht  191  of  oriented  DNA  fibers  by  the  wet  spinning  technique. 
Samples  prepared  in  different  ionic  environments  have  made  it  possible  to  compare 
the  radiolysis  products  of  DNA  obtained  with  different  types  of  radiation  and  with 
different  orientations  of  the  DNA  fiber  axis  relative  to  the  radiation  beam  direction. 
Results  have  been  surprising.  We  note  here  that  exposure  of  such  samples  to  gamma 
radiation  at  77°K  and  measured  at  77°K  resulted  in  EPR  spectra  interpreted  in  terms 
of  G  radical  cations*  and  T  and/or  C  radical  anions  /10,1 1/.  Exposure  of  samples 


S-THYMYL  RADICAL  (TH*)  GUANINE  CATION 

Fig.  1.  Free  radicals  identified  by  EPR  spectroscopy  at  77°K  in  DNA  irradiated 
at  77®K. 

under  similar  conditions  to  neutrons  resulted  in  similar  spectra  if  the  neutron  beam 
was  perpendicular  to  the  DNA  helix  axis;  however,  if  the  beam  was  parallel  to  this 
axis,  the  Em  spectra  indicated  the  presence  of  the  neutral  radiation  product,  the 
dihydrothymidin-5-yl  (TH*)  radical  /12/.  For  proton  irradiation,  however,  the  TH* 


•The  following  symbols  are  used  for  the  DNA  bases:  G=guanine;  A=adenine; 
Ocytosine;  T=thymine 
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spectrum  was  observed  with  both  parallel  and  perpendicular  beam  orientations  /13/. 
Figure  1  summarizes  the  free  radicals  identified  in  irradiated  DNA  at  low  tempera¬ 
tures. 


Attempts  have  been  made  to  explain  these  results  in  terms  of  a  variety  of  microscopic 
long-range  transfer  mechanisms,  such  as  triplet  exciton  migration  /12,14/,  polaron 
transport  /14/  and  solitons  /14/.  We  discuss  these  mechanisms  critically  in  view  of 
recent  research  advances.  We  first  consider  whether  the  occurrence  of  different 
radicals  produced  by  different  kinds  of  ionizing  radiation  are  artifacts  or  properties 
of  the  systems.  We  then  address  whether  the  microscopic  energy  deposition  pattern 
is  sufficient  to  explain  the  results  or  whether  long-range  migration  of  energy  and/or 
charge  is  required  and  how  this  migration  might  best  be  described. 

*  Specifically,  we  question  whether  neutral  triplet  or  singlet  state  migration  provides 
an  explanation  of  the  neutron  results.  This  is  followed  by  a  short  discussion  as  to 
whether  solitons  or  other  DNA  collective  states  or  polaron  transport  are  involved  in 
the  neutron  and  proton  results.  We  then  address  whether  other  microscopic  track 
phenomena  are  operative.  Charge  recombination  and  charge  trapping  in  dense 
particle  tracks  are  briefly  discussed  and  their  effects  on  total  radical  yields  are 
summarized.  The  paper  inconcludes  a  brief  discussion  of  the  reaction(s)  that  might 
account  for  the  appearance  of  TIT  at  low  temperatures  and  concludes  with  sug¬ 
gestions  for  plausible  future  research. 

EXPERIMENTAL  BACKGROUND  SUMMARY 

Samples:  All  experiments  were  performed  at  77°K  with  either  oriented  Na-DNA 
samples  prepared  by  Rupprecht  191  using  the  wet  spinning  technique  or  in  some  of 
the  neutron  experiments,  by  Arroyo  et  al.  / 12/  using  the  procedure  of  Rupprecht. 
Samples  for  the  proton  irradiations  were  carefully  prepared  because  of  the  proton’s 
short-range  in  the  target  medium,  approximately  0.5  mm  for  4  MeV  protons.  For 
experiments  in  which  the  DNA  fibers  were  oriented  perpendicularly  to  the  beam, 
samples  were  made  by  pressing  together  a  sufficient  number  of  sheets  of  oriented 
DNA  to  give  a  sample  thickness  greater  than  the  proton  range.  For  the  parallel 

,  proton  irradiation  configuration,  samples  were  carefully  sliced  from  a  block  of 
oriented  DNA  so  that  the  ends  of  samples  were  not  bent  over  so  as  to  introduce 
spurious  perpendicular  components  into  the  EPR  data.  Although  the  samples  used 

•  for  the  proton  experiments  were  approximately  two  years  old,  they  were  preserved 
under  controlled  conditions.  X-ray  diffraction  studies  on  several  specimens  (personal 
communication,  A.  Rupprecht)  demonstrated  some  distortions  due  to  the  pressure 
required  to  make  the  DNA  film  layers  stick  together  in  the  original  preparation  and 
in  the  original  splicing  of  the  samples,  but  nothing  significant  to  preclude  failure  to 
observed  possible  orientation  effects.  The  DNA  conformation  /15/  as  determined 
from  independent  samples  is  a  mixture  of  the  A-DNA  form,  where  base  planes  are 
at  70°  angles  to  the  helix  axis,  and  the  B-DNA  form  (57%),  where  base  planes  form 
right  angles  to  the  DNA  helix  axis.  All  samples  were  equilibrated  with  water  vapor 
at  75%  relative  humidity;  this  gives  a  moisture  content  sufficient  to  almost  fill  the 
primary  DNA  hydration  shell  /16/.  Figure  2  gives  a  plot  of  the  number  of  water 
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molecules  per  nucleotide  at  •  faction  of  the  relative  humidity  in  the  hydration 
chamber/?/. 

Figure  2  emphasizes  that  even  at  0%  relative  humidity,  there  are  2J  H,0  molecules 
attached  to  die  sodium  phosphate  groups  per  nucleotide.  These  water  molecules  are 
not  removable  via  vacuum  desiccation  of  DNA  samples  /17A  For  our  purpose  we 
note  that  the  tightly  bound  HjO  molecules  (die  first  12*15  water  molecules)  are 
impermeable  to  cations  /18/.  Another  characteristic  of  this  inner  hydration  shell  is 
that  it  does  not  form  an  ice-like  structure.  The  addition  of  approximately  nine  more 
HjO  molecules  per  nucleotide  completes  the  filling  of  DNA  primary  hydration  layer. 
It  is  in  this  latter  layer  of  water  molecules  where  charge  transfer  could  possibly 
occur.  Thus  for  a  relative  humidity  of  75%  the  outer  hydration  layer  is  nearly  filled, 
and  furthermore  contains  sufficient  water  molecules  to  support  charge  transport  along 
the  DNA  backbone. 


Humidity  (%) 


1 0uter 
nyoratton 
She* 


•e 


Inner 


Fig.  2.  Hydration  of  double  stranded  salmon  sperm  DNA  as  a  function  of  relative 
humidity /7/. 
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Radiation  Sources:  The  gamma  radiation  /10.11/  wu  from  a  “Co  source.  The 
neutron  source  was  a  small  TRIGA  research  reactor.  Bismuth  *HielH»ng  was  utilized 
to  reduce  the  gamma  component  to  less  than  3%  of  the  total  dose.  This  procedure 
(unfortunately)  introduces  a  broad  neutron  spectrum  (neutron  with  energies  of  10*eV 
to  more  than  107eV)  with  a  broad  peak  near  1  MeV.  Doses  were  measured  by  an 
indium  foil  technique.  The  principle  neutron  reaction  for  this  energy  range  in  tissue¬ 
like  (water-like)  media  is  elastic  scattering  of  protons.  The  divergence  of  the  neutron 
beam  induces  a  diverging  recoil  proton  beam  suggesting  that  any  orientation  effects 
observed  with  neutrons  might  well  be  enhanced  if  a  well  collimated  proton  beam 
were  employed.  Such  beams  were  obtained  from  a  2  MV  tandem  that 

provided  monoenergetic  protons  with  energy  up  to  4  MeV.  As  protons  have  short 
track  ranges  at  these  energies,  only  the  front  surfaces  of  the  DNA  samples  were 
irradiated.  Table  1  summarizes  the  parameters  for  the  three  experiments. 

Table  1  Irradiation  Parameters 


Investigator 

Radiation 

Dose 

(kGy) 

Dose-Rate 

(kGy/hr) 

Graslund  et 
al.  I4JSJ 

Gamma 

5 

9 

Arroyo  et 
al.  m 

Neutron 

13 

25 

Miller  et 
al.  78/ 

Proton 

30 

130 

EPR  Analysis:  X-band  EPR  spectra  at  100  kHz  and  1  to  2  gauss  modulation 
amplitude  were  obtained  with  microwave  powers  of  <  3pW  (gamma  experiments), 
20  pW  for  spectra  perpendicular  to  the  incident  neutron  beam,  200  pW  for  the 
parallel  radiation  configuration  and  12.3  mW  for  the  proton  experiments.  Lower 
power  was  used  in  the  gamma  experiments  to  avoid  saturation  of  the  ion  radical 
signals. 


RESULTS 

Irradiation  with  gamma  rays  resulted  in  the  spectra  shown  in  Figure  3  /10/.  These 
spectra  were  interpreted  as  indicating  the  presence  of  the  guanine  radical  cation  and 
the  thymine  or  cytosine  radical  anion.  Recent  studies  of  Bernhard  /19.20.21/  indicate 
that  the  anion  probably  is  a  cytosine  and  that  it  may  well  be  protonated  /21/. 
Steenken  Hit  has  stressed  the  importance  of  protonation  in  influencing  electron 
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attachment  For  example,  within  a  base  pair,  proton  transfer  will  stabilize  one 
electron  reduction  of  cytosine.  Bernhard  /19/  has  shown  that  the  probability  of 
electron  attachment  follows  the  relation  T  >  C  »  A  >  O,  which  is  in  agreement  with 
calculated  electron  affinities  of  the  bases  /23,24/.  The  stacking  in  duplex  DNA  in 
aqueous  media,  however,  permits  electron  transfer  between  the  bases,  and  proton 
transfer  horn  a  hydrogen  bonded  partner  may  occur.  Thus  the  tendency  of  the 
hydrogen  bonded  partner  to  deprotonate  becomes  an  important  factor.  The 
observation  that  A  is  a  poorer  proton  doom  than  G  means  that  C"  is  stablized  by 
base  pairing  more  than  T\  Bernhard  /19/  has  shown  that  the  following  relative 
concentration  relationships  are  expected  in  double  stranded  DNA:  (C‘  H*]  »  [A'* 
IT]  -  (Tl  »  [G'l  Note  that  thymine  is  the  only  base  that  should  strongly  form 
an  electron  adduct  but  is  unlikely  to  protonate  reversibly  in  this  environment  On 
heating  the  sample  and  then  recooling  to  77*K  the  EPR  spectrum  shows  the  formation 
of  neutral  TH  radicals  /10,11/,  and  the  disappearance  of  the  ionic  radicals  and  a 
reduction  in  the  number  of  total  spins.  This  is  critically  important  in  the  explanation 
for  the  appearance  of  TH*  radicals. 


Fig.  3.  EPR  spectra  of  calf-thymus  Na  DNA  irradiated  with  gamma  rays  /1(V. 

In  contrast  to  the  gamma-induced  spectra,  proton  irradiation  gave  the  spectra  shown 
in  Figure  4  /13/,  indicative  of  the  presence  of  neutral  TH*  radicals.Spectra  are  similar 
both  for  parallel  and  perpendicular  irradiation  configuration.  For  neutron  irradiation 
the  EPR  spectra  obtained  at  77°K  are  shown  on  Figure  3. 
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Protons  Perpendicular 


Parallel  To  DNA  Fiber  Axis 
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Fig.  4.  EPR  spectra  at  77*K  with  4  MeV  protons  /13/. 

When  the  neutron  beam  was  directed  parallel  to  the  DNA  fiber  axis,  the  neutral 
radical  TH*  was  formed;  however,  when  the  mckfent  neutron  flux  was  perpendicular 
to  the  fiber  axis,  die  spectra  showed  the  presence  of  radical  anions.  The  EPR  data 
raise  die  following  questions:  (1)  by  what  mechanismfs)  does  proton  radiation  give 
the  neutral  protonated  reduction  product  (TIT)  whereas  gamma  radiation  at  77*K 
gives  only  ionic  products?  and  (2)  by  what  mecbanism(s)  does  neutroo  irradiation 
give  spin  products  dependent  on  die  direction  of  die  neutron  beam? 


Are  the  results  with  protons  and  neutrons  due  to  some  problem  with  die  experimental 
procedures?  We  have  conducted  an  extensive  review  of  die  procedures  used  for  the 

neutroo  and  proton  experiments.  Id  the  neutroo  experiments,  two  samples  were 
mounted  side  by  side  within  a  liquid  nitrogen  cooled  inadiatioo  chamber  so  that  the 
fiber  orientation  of  one  sample  was  perpendicular  to  the  neutron  beam  and  the 
orientation  of  the  second  sample  was  parallel  to  the  incident  beam.  The  samples, 
irradiated  simultaneously,  were  treated  identically  so  far  as  is  known.  Furthermore 
die  experimenters  successfully  observed  die  ion  radical  signals  after  gamma 
irradiations,  and  samples  prepared  by  Rupprecbt  and  by  Arroyo  /12/  gave  the  same 
results.  While,  of  course  we  can  not  completely  eliminate  the  possibility  that  errors 
were  made,  we  have  considerable  confidence  in  the  results. 


Fig.  5.  EPR  spectra  from  calf-thymus  Na  DNA  irradiated  with  TRIGA-re actor 
neutrons /12/. 

For  die  proton  study,  individual  samples  were  irradiated  in  a  liquid  nitrogen  cooled 
chamber  and  transferred  to  the  EPR  cavity.  During  the  transfer  the  samples  were 
exposed  to  a  temperature  >  77°K  for  less  than  one  second.  Comparing  this  time  with 
the  times  at  various  temperatures  required  for  the  appearance  of  TH*  /1 1/,  we  feel 
confident  that  this  step  does  not  account  for  the  results.  Furthermore,  the  samples 
were  examined  by  Rupprecht,  and  by  his  criteria,  which  included  determination  of 
X-ray  diffraction  patterns,  they  were  satisfactory  for  die  experiment  Thus  we 
conclude  that  the  results  are  not  a  consequence  of  any  obvious  or  apparent 
experimental  problem. 


In  the  proton  experiments,  are  the  results  due  to  macroscopic  heat  build-up  during 

the  irradiation  accompanied  bv  inefficient  heat  transfer  from  sample  to  coolant? 

Elementary  calculations  show  that  if  all  the  energy  deposited  in  the  sample  were 
retained  and  converted  to  heat,  the  temperature  increase  would  be  of  the  order  of 
10°K.  Using  reasonable  estimates  of  the  cooling,  and  considering  the  dose  rate,  etc., 
we  estimate  a  temperature  difference  of  about  .002°K  between  the  thermal  reservoir 
at  77°  and  the  irradiated  part  of  the  sample.  Such  a  small  temperature  difference, 
even  if  underestimated  by  a  factor  of  103  would  not  be  nearly  enough  to  cause  the 
formation  of  TH*  by  the  reaction  path  observed  by  GrSslund  and  coworkers  /10,1 1/. 


ited  within  the  proton  tracks  in  both  the 
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experiments  explain  the  appearance  of  the  orotonatcd  thymine?  Several  investigators 
/25,26,27/  have  suggested  that  fast  heavy  ions  with  large  stopping  power  result  in 
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transient  localized  bested  regions  around  the  ion  tracks.  To  estimate  this  effect  we 
follow  the  treatment  given  by  Mozutndcr  /28/  wbo  adopted  a  Gaussian  statistical 
approximation.  Solving  the  standard  macroscopic  beat  equation  for  low  LET 
radiation  one  fin ds  the  distribution  of  the  excess  temperature,  where  the  spur  is 
approximated  by  a  sphere  with  size  parameter  r„  to  be: 

AT(1)(r,t)  =  T0  (1  +  4  WrtY*2  expt-i2/^  +  48t)} 


Here  T.  is  the  maximum  excess  temperature  at  the  spur’s  center,  r  is  the  distance 
from  the  spur  center,  and  8  is  the  thermal  diffusivity  and  equals  X/p the  heat 
conductivity  (X)  divided  by  the  product  of  the  medium  density  (p)  and  specific  heat 
(C,)  at  constant  volume.  T,  is  determined  by  the  relationship: 

oo 

E  =  /  pCv  T0  4jc  r2  exp(-r2/r^)dr 
Jo 

where  E  denotes  the  average  energy  deposited  per  charge  pair  in  the  spur,  -30eV. 
For  high  LET  ionizing  particles  traversing  s  medium  s  large  number  of  phonons  are 
created  resulting  in  a  high  density  of  transient  temperature  pulses.  Ihese  heated 
regions  overlap  and  thereby  form  a  cylindrical  excess  thermal  distribution  with  initial 
radius  r,  that  can  be  approximated  as: 

AT^r.t)  =  T0(1  +  45t/r|)-1  exp{-r7(r§  +  48t)} 


where  r  is  die  transverse  distance  from  the  trade  axis.  If  S  denotes  the  energy  loss 
per  unit  distance,  then 

OO 

s  =af  =  f  pcv  To2jcr  expt-^/r^dr 

determines  the  initial  temperature  T,  at  the  center  of  the  cylindrical  trade.  Figures 
6A  and  (SB  illustrate  the  tune  dependence  of  local  heating  for  several  r  values 
assuming  8  «  10rScm’/sec,  r.  *  20A,  E  *  30eV  and  S  *  5eV/A. 

Variations  for  aT  are  quite  similar  and  in  both  cases  decay  to  the  ambient 
temperature  in  approximately  Iff10  to  10*  seconds.  Although  these  calculated 
temperature  transients  are  similar,  the  effect  of  the  heat  spike  on  reaction  rates  is 
strongly  dependent  on  dimensionality.  Assuming  an  activation  energy  of  8 
Kcal/mole,  a  substrate  concentration  of  10“  molecules/cm*  and  a  collision  frequency 
of  lO"  s*1 
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Fig.  6.  Excess  temperature  transients  (A)  AT(,)(r,t)  in  a  spur  Qow-LET  radiation) 
and  (B)  AT^r.t)  in  a  cylindrical  track  (high-LET  radiation)  at  several  distances 

( - ,  r  *  lOA, - ,  r  *  15  A,  and - ,  r  *  20  A)  from  the  center  of  the 

spur  or  the  cylindrical  track  axis. 

K  =  10  "expf-SOOO/RT)  cmWmolecule 


gives  a  low  LET  reaction  time  of  *  0  Jpsec  (assuming  T,  •  400°K).  For  high  LET 
izradiatioa  with  S  *  SOOeV/A  corresponding  to  a  T,  of  104  °K  (RT*  0.83eV)  the 
reaction  time  is  estimated  to  be  1 5  x  Wsec,  a  time  comparable  to  the  duration  of 
the  thermal  pulse.  This  simple  calculation  demonstrates  that,  at  least  for  low  LET 
radiation,  the  induced  transient  increase  in  the  local  temperature  does  not  effect 
chemical  processes.  This,  however,  need  not  be  valid  for  very  high  LET  radiation 
since  the  temperature  pulse  duration  can  be  comparable  to  chemical  reaction  times. 

Similar  calculations  for  4  MeV  proton  tracks  yields  an  initial  temperature  of  about 
100°K,  which  is  probably  not  sufficient  to  stimulate  the  formation  of  TH*. 
However,  errors  in  the  estimates  could  well  bring  the  temperature  into  the  necessary 
range  of  about  200*K.  Thus  intra-track  heating  by  proton  radiation  during  and 
shortly  thereafter  can  not  be  excluded  as  a  possible  mechanism  involved  in  the 
formation  of  TH*. 
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Several  microscopic  processes  have  been  suggested  to  explain  the  observed 
anisotropy  in  the  neutron  experiments.  We  discuss  briefly  whether  any  of  the  quasi¬ 
particles,  excitons,  solitons  or  polarons  or  other  track  phenomenon  can  account  for 
the  anisotropic  neutron  EPR  spectra. 

Does  the  migration  of  triplet  and/or  singlet  excitons  provide  an  explanation  of  the 
neutron  EPR  spectra?  In  the  original  paper  by  Arroyo  et  al  /12/  it  was  suggested  that 
the  migration  of  triplet  states  might  provide  an  explanation  of  the  neutron  EPR  data. 
We  now  believe  this  transport  mode  cannot  provide  the  large  asymmetry  in  the 
transport  of  energy  in  DNA  (see  the  discussion  on  microscopic  intra-track  spike 
model)  needed  to  account  for  the  neutron  data.  Both  singlet  and  triplet  exciton 
migration  (within  oriented  DNA)  is  probably  limited  to  at  most  10  base  pairs  due  to 
their  short  life  time  and  extensive  trapping  expected  to  exist  in  highly  disordered 
systems.  Although  excess  hole  and  electron  band  widths  along  t  le  DNA  sugar- 
phosphate  backbone  are  estimated  to  “400  cm'1  and  1200  cm'1,  respectively  1291,  the 
neutral  excited  state  is  considerably  narrower,  about  10  to  70  cm'1  1301.  These 
theoretical  values,  although  admittedly  crude,  support  the  contention  that  mobile 
neutral  excited  states  have  very  short  transfer  distances.  Numerous  experiments 
support  this  viewpoint.  As  an  example  we  note  that  exciplex  formation  in  DNA  is 
a  rapid  process  /3 1/  estimated  to  be  the  order  of  10l2sec'';  thus  excited  singlet  transfer 
needs  to  occur  in  less  than  10',2sec.  Denoting  the  mean  transfer  distance  by  d,  we 
have 


d82  =  2Fa2tT 


if  only  nearest  neighbor  transfers  are  considered  with  a  Fdrster  rate  F  «  1013  sec1. 
Here  a  is  the  neighboring  base  pair  stacking  distance  (G.34nm  in  B-DNA)  and  tx  is 
the  trapping  time.  Assuming  an  exciplex  formation  time  of  1012  sec  gives  a  transfer 
distance  d,  of  1.5  nm  or  approximately  4  to  5  base  pairs.  Triplet  state  transfer  has 
been  investigated  by  Eisinger  and  Lamola(32)  and  by  Gueron  and  Shulman  (33)  to 
note  but  a  few.  Probably  the  most  definitive  experimental  evidence  for  triplet 
migration  in  DNA  has  been  provided  by  measuring  the  quenching  of  DNA  phospho¬ 
rescence  by  metal  ions.  Analysis  of  the  data  by  Isenberg  and  coworkers  (34)  on  Ni2\ 
Co2+  and  Mn2+  quenching  gives  triplet  transfer  distances  comparable  to  that  of  the 
singlet  state,  (0.8  to  1.5  nm).  We  therefore  can  safely  conclude  that  neutral  state 
migration  range  is  short  in  DNA  and  does  not  constitute  an  effective  mechanism. 

Does  a  microscopic  temperature  spike  model  using  theoretical  track  codes  account 
for  the  anisotropy  in  neutron  EPR  data?  The  most  obvious  difference  in  the  pattern 
of  energy  absorption  in  the  two  irradiation  geometries  is  that  charged  particles 
traversing  the  sample  nearly  parallel  to  the  helical  DNA  axis  have  a  greater 
probabliltiy  for  multiple  energy  transfers  to  the  same  DNA  molecule  than  do  particles 
incident  on  the  sample  perpendicular  to  the  fiber  orientation.  This  difference  in  the 
spatial  pattern  of  energy  absorption  should  have  no  effect  on  the  production  of  free 
radicals  if  energy  and/or  charge  move  between  different  DNA  fibers  as  freely  as  they 
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are  transported  along  a  single  DNA  chain.  The  macroscopic  model  developed  by 
Miller  et  al.  /35/  (the  quasi  particle  being  transported  is  not  identified)  assumes  that 
an  orientation  dependence  of  radical  yields  is  indicative  of  intramolecular  energy 
and/or  charge  transfer  in  DNA.  The  asymmetric  3-dimensional  heat  equation  is 
solved  in  the  limit  where  only  the  shortest  transverse  relaxation  time  (t,  *  b2/5.8D„ 
where  b  is  the  DNA  radius  and  D,  is  the  transverse  diffusioL  coefficient)  is  included. 
Miller  et  al.  showed  that  in  this  approximation  the  local  temperature  at  energy 
deposition  site  x*  (determined  by  Monte  Carlo  track  simulation)  and  at  time  t  is  given 
by 


T(x,t)  =  Ts  +  exp(-t/x,  )Xjk(x,t/T|) 

and 


Tk(x,t)  =  Tk0  (1  +  t/T|)  1/2exp{~2^1  ^ 

where  T,  is  the  ambient  sample  temperature,  x,  =  A/2D,  is  the  longitudinal  thermal 
relaxation  time,  D,  is  the  longitudinal  diffusion  coefficient,  A  is  the  half-width  of  the 
superexcited  state  (taken  to  be  3.4A,  the  nearest  neighbor  base-pair  distance)  and 

Tko  =  £fc/b2ApC(27C3)1/2 

where  C  and  p  are  the  specific  heat  and  sample  density  and  ek  is  the  absorbed 
energy  at  site  x*. 

An  essential  difference  between  the  Henriksen  et  al  /26,27/  application  of  Mheat 
waves"  to  calculate  conversion  of  the  primary  radicals  in  proteins  and  the  thermal 
spike  model  of  Miller  and  coworkers  /35/  is  that  of  reformulating  the  thermal  spike 
model  in  terms  of  track  structure  rather  than  simply  stopping  power.  In  view  of 
recent  findings  by  Bernhard  /1 9,20,21/,  electrons  liberated  in  the  decay  of  superex¬ 
cited  stales  by  autoionization  (see  Figure  7)  eventually  form  C'. 

Energy  from  other  deposition  events  in  the  same  DNA  chain  are  treated  as 
temperature  spikes  that  spread  throughout  the  sample  asymmetrically.  Thermal 
diffusion  along  DNA  chains  is  assumed  to  be  1000  fold  higher  than  the  transverse 
thermal  diffusion.  Although  at  present  we  are  unware  of  any  experimental  evidence 
to  support  this  large  value,  it  should  be  noted  that  anisotropy  in  electrical  conduc¬ 
tivity  greater  than  100  has  been  reported  for  organic  systems  such  as  polydiacetylene 
/36/.  The  longitudinal  diffusion  coefficient  was  taken  to  be  comparable  to  H30  ( 10' 3 
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Fig.  7.  Possible  decay  modes  of  energy  absorbed  from  ionizing  radiation  in 
hydrated  DNA  samples. 

cm2  sec*1)  and  protonation  of  the  thymine  radical  was  approximated  by  quasi  first- 
order  kinetics  with  an  activation  energy  between  0.2  and  0.6  eV.  The  mechanism  for 
conversion  of  primary  radical  anions  to  TH\  which  probably  involves  electron 
transfer  from  C  to  T  followed  by  irreversible  protonation  at  C6,  is  assumed  to  be  the 
same  as  that  which  operates  in  thermal  annealing  experiments  /1 1/.  Although  the 
model  does  predict  greater  conversion  of  primary  radical  anions  to  TH-  for  a  proton 
flux  oriented  parallel  to  the  DNA  fibers,  it  cannot  account  for  the  3  fold  greater  total 
radical  yields  observed  in  the  parallel  neutron  irradiation  configuration  without 
additional  considerations  of  electron  trapping  and  ion  recombination.  Furthermore, 
the  model  predicts  that  both  primary  radical  anions  and  TH*  should  be  seen  in  the 
parallel  case,  which  does  not  agree  with  the  observations  reported  by  Arroyo  et  al. 
/12/  where  radical  anion  species  were  detected  in  the  only  in  the  perpendicular  case. 

Are  solitons  involved  in  energy  transfer  in  DNA?  Solitons  have  been  proposed  as 
a  mechanims  of  long-range  energy  transfer  in  DNA  /14, 37,38/.  These  collective 
states  have  been  invoked  to  explain  the  hydrogen-deuterium  exchange  reaction 
observed  in  double-stranded  polynucleotides  /39,40/;  however,  problems  with  this 
explanation  have  been  noted  by  Benight  et  al.  /41/.  There  is  currently  no  generally 
accepted  theory  of  solitary  waves  in  DNA.  As  a  first  approximation,  one  might 
visualize  a  DNA  solitary  wave  as  a  configuration  formed  by  partial  underwinding  of 
the  helix  over  several  base  pairs  with  some  disruption  of  interstrand  hydrogen  bonds. 
Yomosa  1421  developed  a  simple  DNA  soliton  model  that  allows  one  to  estimate  their 
lowest  energy  Eq,  size  Lq,  and  velocity  v0  by  fitting  model  Hamiltonian  calculations 
to  the  temperature  dependence  of  the  equilibrium  constant  for  hydrogen-deuterium 
exchange  / 40 /.  Although  details  of  the  Yomosa  model  of  open  states  are  likely  to 
be  incorrect  (e.g.  the  assumption  of  complete  base-pair  rotation  about  the  phosphate- 
sugar  backbone  and  the  neglect  of  DNA  backbone  bending  energy),  it  does  provide 
reasonable  estimates  of  the  above  parameters.  Under  the  assumption  that  base 
stacking  is  preserved  in  soliton  excitation  and  with  a  Hamiltonian  which  includes 
interstrand  hydrogen  bonding,  intrastrand  stacking,  and  torsional  energy,  Yosoma  742/ 
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calculated  Eg  *  0.35  eV,  which  is  sufficient  to  overcome  the  activation  barrier  to 
irreversible  protonation  of  pyrimidine  radical  anions  at  C6.  Lg  =  10  base  pairs,  which 
demonstrates  the  non-local  character  of  this  excitation  mode,  and  v0  *  8.3  x  10s 
cm/sec,  which  is  less  than  the  velocity  of  sound  in  DNA  as  it  should  be. 

The  soliton  migration  distance  d  can  be  estimated  by  noting  that  in  the  radiation- 
induced  thermal  spike  model  of  Miller  et  al.  /35/,  the  heat  wave  has  a  lifetime  of 
approximately  1  ns,  a  value  considerably  shorter  than  the  measured  open-state 
lifetime  of  10"7  s  /43/;  thus  d  *  83  nm.  Hence,  intramolecular  energy  transfer  by 
solitons  is  a  credible  mechanism  for  the  necessary  long-range  energy  transfer  inferred 
from  Monte  Carlo  calculations  /35/  to  explain  the  neutron  data.  Unfortunately,  a 
soliton  model  without  charge  recombination  and  trapping  cannot  account  for  the 
observed  neutron-induced  radical  yields.  This  limitation  might  be  removed  by 
inclusion  of  the  effects  of  soliton  pinning;  however,  a  theory  of  this  process  in  DNA 
is  not  currently  available. 

Is  polaron  transport  involved  in  the  neutron  or  proton  results?  It  is  possible  that 
long-range  migration  (LRM),  either  by  electrons  or  protons  could  account  for  the 
neutron  data.  Charge  transfer  in  organic  systems  is  well  documented  /30/.  Transport 
through  stacked  DNA  bases  has  been  shown  by  van  Lith  and  coworkers  /44/  to  be 
the  order  of  9  nm,  which  is  not  long  enough  to  explain  the  presence  of  TH*  radicals 
at  77°K.  Futhermore,  these  investigations  demonstrated  that  transport  through  the 
phosphate-sugar  backbone  is  also  not  sufficient  Experimental  data  do,  however, 
support  the  possibility  of  LRM  within  the  structured  water  layers  of  DNA  hydration. 
Using  transient  microwave  techniques,  van  Lith  et  al.  /44/  observed  a  mobility  p  of 
2. 5x10°  m2  V'1  sec'1,  which  they  associate  with  transfer  of  "dry"  electrons  in  the 
water  layer,  provided  the  water  concentration  was  above  the  critical  weight  fraction 
F0  *  0.44.  For  comparison,  mobilities  of  2x10*  m2  V1  sec'1  and  lxlO4  m2  V'1  sec'1 
have  been  reported  for  horizonally  stacked  phthalocyanine  columns  /45/  and 
polyethylene  /46/,  repectively.  Above  the  critical  water  content,  the  conductivity  is 
linear  in  F.  If  the  carrier  lifetime  (determined  primarily  by  trapping  at  defects)  is  an 
order  of  magnitude  longer  than  the  unsolvated  electron  in  pure  ice,  then  the  mean 
transfer  distance  at  -78°  C  is  about  100  nm.  More  recent  estimates  by  Warman  et 
al.  /4 5/  give  a  maximum  transfer  distance  of  only  115  base  pairs  or  about  40  nm. 

Although  electron  transfer  is  sufficient  to  explain  LRM  in  DNA,  proton  transport  is 
another  posibility.  Kunst  and  Warman  /47/  attributed  the  longer-lived  conductivity 
transient  for  ice  (after  electrons  are  trapped)  to  mobile  protons.  If  this  is  the  case, 
then  narrow-band  polaron  transport  /48 /  is  the  appropriate  theoretical  formalism  to 
describe  the  dependence  of  mobility  on  temperature.  In  this  theory 

p  =  (ea2/kT)  x1/2v0  F(T)1/2exp(-F(T)) 


with 
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F(T)  =  2Wp/hv0  csch(hv</2kT) 

where  w,  is  the  poUroo  binding  energy,  a  is  the  interatomic  distance,  and  ve  is  the 
vibrational  frequency  of  the  phonon  that  interacts  most  strongly  with  the  carrier. 
Using  these  equations,  we  estimate  that  the  proton  mobility  in  hydration  layers  of 
DNA  is  between  0.1  and  1  cm2  V'1  sec'1  at  77°K.  To  achieve  a  migration  distance 
of  lOOnm,  protons  with  mobilities  in  this  range  would  require  lifetimes  of  4  to  40  ns, 
which  are  consistent  with  existing  results  on  the  trapping  of  electrons  and  protons  in 
ice  /47/.  Hence,  mechanisms  of  charge  transport  exist  wife  sufficient  range  to 
account  for  fee  orientation  effect  on  radical  types  reported  for  neutron  irradiation  but, 
as  was  the  case  for  the  thermal  spike  model,  trapping  and  recombination  must  be 
included  in  any  model  of  radical  yields. 
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data?  As  noted  above,  considerations  of  trapping  and  recombination  in  addition  to 
LRM  of  any  quasi-particle  are  needed  to  explain  the  factor  of  3  difference  between 
total  radical  yields  observed  in  parallel  and  perpendicular  neutron  irradiation  l\2J. 
Miller  and  Swenbeig  738/  have  suggested  ways  to  model  this  effect  as  well  as  the 
absence  of  TH-  radicals  in  the  perpendicular  configuration.  The  observations  of  van 
Life  et  al.  /44/  of  electron  LRM  in  hydrated  DNA  following  nanosecond  pulses  of  3 
MeV  electrons  suggest  that  autoionization  of  superexcited  states  (as  illustrated  in  Fig. 
7)  produces  quasi-free  electrons  in  the  hydration  layers  of  DNA  where  they  undergo 
LRM  before  relaxation  into  another  trap.  Films  of  oriented  DNA  may  be  similar  in 
some  repects  to  the  qua&i-one-dimensional  semiconductors  illustrated  schematically 
in  Figure  8,  below.  In  the  perpendicular  irradiation  configuration,  yields  of  primary 
radical  anions  and  cations  are  mainly  determined  by  the  competition  between 
geminate  recombination  and  electron  trapping  at  preexisting  defects  (denoted  by 
squares  in  Fig.  8).  In  the  parallel  case,  the  high  mobility  path  of  electrons  (for 
polydiacetylene  mobility  parallel  to  the  fiber  direction  can  be  100  times  greater  than 
mobility  perpendicular  to  the  fibers  /3d/)  contains  trapped  positive  ions  that  will 
reduce  the  radical  yield  by  nongeminate  recombination  of  the  ionic  precursors. 
Solution  of  the  coupled  differential  equations  that  describe  radical  production  and 
decay  in  this  model  will  be  given  in  a  future  publication.  Here  we  note  only  that, 
in  contrast  to  the  solitoo  model  depicted  in  Figure  7  where  collective  vibrational 
excitations  are  fee  mobile  low-lying  states,  ejected  electrons  are  fee  quasi-particle 
wife  asymmetic  mobility  in  the  model  depicted  in  Figure  8.  Vibrational  excitations 
that  accompany  both  the  production  and  decay  of  ion  pairs  are  assumed  to  be 
localized. 


Our  analysis  indicates  feat  no  energy  or  charge  transport  model  can  fully  account  for 
anisotropy  (in  the  neutron  experiments)  observed  in  the  EPR  data  if  the  effects  of 
carrier  trapping  and  recombination  are  neglected.  Geminate  recombination, 
applicable  in  the  low  density  limit  where  an  isolated  pair  of  oppositely  charged 
geminate  particles  can  be  considered  as  an  isolated  system,  was  developed  in  the  late 
1930’s  by  Onsager  749,50/  and  was  designed  as  the  steady  state  solution  of  the 


(10)196 


C  B.  Swuti  tt  d. 


ONAnSCM  OMCNTATKM 


Fig.  8.  Schematic  diagram  of  electron  trapping  and  recombination  in  a  quasi  one 
dimensional  semiconductor.  Squares  and  circles  denote  preexisting  and  radiation- 
induced  traps,  respectively.  Broken  lines  represent  electron  trajectories.  (Redrawn 
from  /36/). 

charge  pair  moving  in  a  condensed  phase  continuum  in  the  presence  of  an  external 
electric  field.  This  theory  has  been  widely  and  successfully  applied  to  photogenera¬ 
tion  of  carriers  in  organic  materials  (see  r*ope  and  Swenberg  /30,51/  for  a  review). 
Within  the  past  decade  the  theory  has  been  extended  by  Hong  and  Noolandi  /SI/  to 
include  transient  effects  with  the  refinement  of  replacing  the  original  Onsager 
assumption  of  a  point  sink  at  the  origin  by  a  recombination  sphere  of  finite  radius 
and  recombination  velocity.  In  addition,  the  seminal  theory  of  Scher  and  Rackovsky 
/52 /  in  which  the  effects  of  the  lattice  and  the  microscopic  molecular  process  are 
considered  has  demonstrated  the  important  role  of  competing  processes,  such  as  the 
decay  rate  of  the  electronically  excited  precursors  to  the  ion-pair  state  that 
dissociates.  This  improved  theory  shows  how  the  two  parameters  of  Onsager’s 
theory,  the  electron  thermalization  distance  and  the  initial  quantum  yield  of  geminate 
pair  generation  can  be  evaluated  if  an  assumption  is  made  regarding  the  initial 
distribution  of  geminate  pair  distances.  Recently,  within  the  context  of  Onsager’s 
continuum  theory,  Mozumder  15 41  has  illustrated  by  extensive  numerical  calculations 
how  the  fractal  geometry  of  the  lattice  can  influence  the  geminate  escape  probability, 
mean  recombination  time,  and  the  reaction  rate.  All  these  theoretical  treatments  are 
applicable  only  in  the  low  ionization  density  limit;  in  the  cases  where  high  LET 
particles  are  involved  the  free  electrons,  ions,  and  radicals  in  the  track  are  so  high  in 
density  that  non-geminate  processes  are  dominant  (see  Schott  /55/).  When  ionization 
is  dense,  no  similar  results  as  found  for  geminate  processes  are  known.  This  is  due 
to  two  complicating  factors;  the  repulsion  between  electrons  and  the  inherent  problem 
of  tracking  the  amount  of  cation  charge  remaining  as  electron  and  cation  are 
neutralized.  Both  of  the  effects  have  been  treated  in  a  mean  field  approximation 
although  one  can  easily  convince  oneself  that  neither  a  mean  field  approximation  nor 
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a  perturbation  approach  can  be  successfully  employed.  Sano  and  Baird  1561  have 
considered  the  model  system  of  two  electrons  in  the  field  of  a  divalent  cation  (  a 
model  analog  of  the  helium  atom  in  atomic  theory).  Even  here  the  results  are 
mathematically  formidable  and  require  numerous  numerical  calculations.  It  seems 
that  the  only  viable  approach  is  to  employ  a  variational  approach,  and  this  obviously 
requires  good  theoretical  insight  on  the  form  for  the  multidimensional  distribution 
functions. 


CONCLUSIONS 

Several  mechanism  of  energy  or  charge  transport  in  hydrated  DNA  have  been 
proposed  that  have  sufficient  range  to  couple  energy  deposition  events  in  DNA  chains 
in  ways  that  may  account  for  the  orientation  dependence  of  radical  yields  in  oriented 
DNA  exposed  to  neutrons  /12 /.  The  details  of  these  mechanisms  will  be  dependent 
on  the  identity  of  the  primary  radiation-induced  species;  however,  the  concepts  of 
relating  orientation  effects  to  energy  or  charge  transport  is  independent  of  whether 
thymine  or  cytosine  is  the  predominant  type  of  electron  gain  center.  The  pattern  of 
energy  deposition  events  in  an  oriented  DNA  chain  interacting  with  protons  depends 
upon  the  velocity  of  the  proton  (magnitude  and  direction)  relative  to  the  helical  axis 
but  should  be  the  same  for  a  primary  or  secondary  flux.  Hence,  the  inconsistancy 
between  experiments  with  neutrons  /12/  and  direct  proton-beam  irradiation  /13/  is 
currently  the  main  impediment  to  understanding  the  unusual  results  observed  with 
neutrons.  Determining  the  reproducibility  of  these  experiments  is  clearly  the  most 
important  next  step;  however,  a  detailed  analysis  of  the  spatial  distribution  of  energy 
deposited  in  DNA  chains  exposed  to  neutrons  from  the  TRIGA  reactor  would  also 
be  helpful.  If  the  dispersion  of  secondary  protons  makes  the  patterns  of  energy 
deposition  in  DNA  chains  essentially  independent  of  their  orientation  relative  to  the 
neutron  flux,  then  models  for  the  orientation  effects  reported  by  Arroyo  et  al.  /12J 
that  are  based  on  an  assumed  difference  in  this  pattern  are  obviously  inappropriate. 
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